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DROUGHT RESPONSE OF TWO MEXICAN OAK SPECIES, 
QUERCUS LACEYI AND Q. SIDEROXYLA (FAGACEAE), IN RELATION 
TO ELEVATIONAL POSITION! 
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To investigate the relationship between the altitudinal distribution of Quercus laceyi and Q. sideroxyla and their physiological 
responses to drought, we measured relative water content (RWC), water potentials (Wpredawn and ‘Pmidday), photosynthesis (Amax), 
stomatal conductance (g), chlorophyll fluorescence (F\/F'm), and spectral reflectance (400-1100 nm) five times during a 7 wk 
acute drought. Quercus laceyi was drought tolerant, while Q. sideroxyla was a drought avoider; Q. laceyi tolerated lower RWC 
(Q. sideroxyla = 54%, Q. laceyi = 44%), Ppa (Q. sideroxyla =—2.6 MPa, Q. laceyi=—3.3 MPa), and ma (Q. sideroxyla = —4.5 
MPa, Q. laceyi=—6.6 MPa). The Fy/Fm also declined first in Q. sideroxyla in wk 6, whereas Fy/F'm did not decline in Q. laceyi 
until wk 7. Amax and g fell in wk 4, 6, and 7 in drought seedlings of both species, suggesting a decline in CO2 assimilation during 
the drought. Leaf spectral reflectance increased with time in response to decreases in leaf photosynthetic pigment concentrations 
in latter weeks of the drought. The results suggest a close association between the altitudinal distributions of these species and 


their adaptation to water stress. 
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The genus Quercus (oaks; Fagaceae) is one of the most 
diverse groups of temperate trees, with greater than 500 species 
in the temperate and subtropical latitudes of the Old and New 
worlds (Nixon, 1993). Mexico represents the origin of oak 
diversification in North America and contains the highest oak 
diversity in the world (Muller and McVaugh, 1972; Rzedow- 
ski, 1978; Gonzalez, 1993; Nixon, 1993). Oak and pine—oak 
forests cover almost half (45%) of the temperate coniferous 
forest of Mexico (World Forest Institute, 1994) and are 
globally significant centers of vascular plant diversity and of 
endemism in both plant and animal species (Manuel-Toledo 
and Jestis-Orddéfiez, 1993). However, the mechanisms under- 
lying oak diversification across Mexico are poorly understood. 

With 63% of its surface area covered by desert, Mexico 
contains the most arid regions in North America. In these 
extreme environments, water is the predominant factor limiting 
carbon gain and primary production (Noy-Mier, 1973; DeLucia 
and Schlesinger, 1991; Lajtha and Getz, 1993; Kozlowski and 
Pallardy, 1997). Yet, no published information quantifies the 
variability of oak drought tolerance in relation to species 
distribution patterns in Mexico. 

Oaks are some of the most drought tolerant trees of the 
temperate latitudes (Abrams, 1990; Ogaya and Pefiuelas, 2003; 
Aranda et al., 2005). North American oaks are adapted to 
drought prone sites through a variety of mechanisms including 
drought avoidance and drought tolerance (Abrams, 1990, 
1996), and the distribution of oak species in the United States is 
tightly linked to the physiological, anatomical, and morpho- 
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logical variation among species and across sites (Kubiske and 
Abrams, 1992; Ashton and Berlyn, 1994). 

Understanding the relationship between oak physiology and 
moisture regimes is especially important in arid regions of 
North America, where precipitation is already greatly exceeded 
by evaporative demand and where the water deficit might be 
increased as a result of global warming (Houghton et al., 2001; 
Loik et al., 2004). Some reports describe the response of 
several oak species in the Rocky Mountains (Neilson and 
Wullstein, 1987; Dickson and Tomlinson, 1996; Ehleringer 
and Phillips, 1996; Balok and St. Hilaire, 2002), yet to our 
knowledge, no prior research explains the effects of moisture 
regime variability across mountain ranges on oak ecophysiol- 
ogy and adaptations to water stress in the oak species 
complexes in Mexico. 

The objective of this study was to measure the physiological 
responses of a high and a low elevation endemic Mexican oak 
species of the genus Quercus to an experimental acute drought 
similar to that experienced by seedlings in the field before the 
onset of the North American Monsoon System (NAMS), which 
occurs in the spring across the southwestern United States and 
Northern Mexico (Adams and Comrie, 1997). Our intent was 
to investigate the physiological mechanisms responsible for the 
distribution of two oak species, Quercus laceyi Small 
(subgenus Leucobalanus; Fagaceae) and Q. sideroxyla Humb. 
et Bonpl (subgenus Erythrobalanus; Fagaceae), on different 
sites within the same mountain range. 

Seedling performance was examined under two moisture 
regimes (unlimited water and drought) on two species 
inhabiting sites that vary dramatically in local environmental 
conditions. Quercus laceyi is the lowest elevation tree of any 
species in the Sierra del Carmen, Coahuila, Mexico, and it is 
the dominant low-elevation oak in savannahs and riparian 
areas. It is a medium stature tree that ranges from 3-8 m in 
height, and it typically colonizes hot, dry sites. Quercus 
sideroxyla is the dominant high-elevation oak in the Sierra del 
Carmen. This species colonizes cooler, wetter sites compared 
to Q. laceyi. It coexists with a mixture of other oak and conifer 
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Fig. 1. 
Mexico. 


Location of the Maderas del Carmen study area in Coahuila, 


species in closed canopy forests that reach over 30 m in height. 
The two species have overlapping elevational ranges, but occur 
in different forest types giving them spatially segregated 
distributions. 

By measuring a suite of physiological parameters during an 
imposed 7-wk drought, we tested the hypothesis that these 
species have distinct physiological responses to drought that 
reflect their different altitudinal distributions. Differences in the 
physiological response of these species to drought may explain 
the niche differentiation or specialization of oaks on different 
sites within the same mountain range (May and MacArthur, 
1972; Pianka, 1976). Furthermore, this research has the 
potential to elucidate the mechanisms responsible for the high 
oak diversity in northern Mexico. 


MATERIALS AND METHODS 


Experimental conditions and experimental design—Ripe acorns were 
collected from 20 mature trees of each species in the Sierra del Carmen of 
Coahuila, Mexico in November 2004 (Fig. 1). In montane woodlands of this 
region, acorns mature in the fall, germinate, and establish under frequent winter 
precipitation, and are then exposed to a prolonged drought from March until 
mid-July when precipitation returns with the NAMS. While stratification is 
necessary for the germination of some oaks of the subsection Erythrobalanus, 
oaks in the southwestern United States and Mexico germinate readily after 
maturation (Nyandiga and McPherson, 1992). 

Six hundred acorns were submerged in water and placed in a 5°C 
refrigerator for 24 h to promote germination. Acorns were planted in a 
Metromix seedbed (Scott’s Co., Marysville, Ohio, USA) inside a greenhouse 
with a daily temperature range of 18-42°C and mean relative humidity of 20% 
(12-45% range) in New Haven, Connecticut on 22 March 2005 (lat. 41°16’ N; 
long. 72°55’ W). Air temperature and relative humidity were measured with a 
Hobo Micro Station Data Logger (Onset Computer, Bourne, Massachusetts, 
USA). Photosynthetically active radiation (PAR) was measured at 30-min 
intervals across the greenhouse bench used in the experiment with a HOBO 
Photosynthetic Light (PAR) Smart Sensor (Onset Computer) attached to the 
data logger. The mid-day peak in PAR varied from 1100 to 1561 


Eighty percent of the acorns planted in our study germinated under 
greenhouse conditions. No information exists that describes the germination 
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success of these two species in the field. However, two similar species, Q. 
arizonica Sarg. (subgenus Leucobalanus) and Q. emoryi Torr. (subgenus 
Erythrobalanus), also in the Sierra del Carmen, had 72% and 21% germination, 
respectively, under field conditions in Arizona (Nyandiga and McPherson, 
1992). 

After 3 wk, when the radical was 3 cm long, germinants were transplanted to 
10-cm plastic pots. Seedlings were transplanted again to 18-cm plastic pots 
filled with Metromix growing medium on 18 September 2005. A modified 
Hoagland’s nutrient solution (Hoagland and Arnon, 1950) was applied monthly 
for 3 mo to all plants before the onset of the experiment. Metromix growing 
medium was chosen because it is a commercially available medium that drains 
well. The soils of the natural habitats of the two species are generally shallow to 
moderately deep and are volcanic in origin. The underlying rocks are 
predominantly extrusives, consisting of lavas and pyroclastics. 

On 21 March 2006, 120 healthy seedlings of each species (240 total) were 
divided into two treatments, drought and control, using a completely 
randomized sampling design. The date of experiment onset was chosen to 
simulate the drought experienced by these species in the field before the onset 
of the NAMS (Adams and Comrie, 1997). High temperatures and _ little 
precipitation characterize the natural drought. No climatic data exist for the 
Sierra del Carmen, but mean annual precipitation for the Chisos Basin 
approximately 150 km to the north is 70 cm (range 32-135 cm from 1948— 
2003). Mean January precipitation is 1.5 cm (range 0-2.5 cm), mean 
precipitation in March before the NAMS is 0.1 cm (range 0-3.7 cm), and 
mean precipitation in July during the NAMS is 8.0 cm (range 0.2—20.5 cm). 
Mean monthly minimum temperatures are 1.8°C in January, 6.3°C in March 
and 17.0°C in July. Maximum temperatures are 14.1°C in January, 20.9°C in 
March, and 29.1°C in July. Relative humidity during the drought (March—July) 
fluctuates between 10 and 25%. Daily greenhouse temperatures and relative 
humidity were not controlled during the experiment, although the mean relative 
humidity and temperatures were within the ranges experienced by seedlings in 
the field during the natural drought. 

The experimental unit was a single seedling in a pot, and soil moisture was 
monitored using gravimetric methods. All pots were watered to saturation (field 
capacity) at the onset of the experiment and were weighed. Control seedling 
pots were watered to 100% of field capacity three times per wk by bringing 
them to their field capacity weight. After initial watering, seedlings in the 
drought treatment were exposed to an acute drought, where water was 
completely withheld for the duration of the experiment. Drought seedlings were 
weighed weekly to determine their percentage of field capacity mass. Drought 
seedlings reached 10% field capacity at the end of the experiment (reached 48 d 
after the onset of the drought treatment). 


Water relations—Predawn (0400 to 0500 hours) and midday (1200 to 1300 
hours) water potentials (pa and Pima) (MPa) were measured at the onset of the 
experiment, and at 2, 4, 6, and 7 wk into the experiment. Three individuals of 
each species-treatment combination (NV = 12) were destructively sampled at 
predawn for the determination of ‘Ypq using a pressure chamber (PMS 
Instruments, Corvallis, Oregon, USA). Water potential measurements were 
repeated at 1200 hours on the same days on 12 additional individuals to obtain 
midday ‘Wma values. 

Leaf relative water content (RWC) was determined for five leaf discs from 
leaves on the same seedlings used in the ‘¥ measurements. Leaf discs were cut 
and immediately weighed on a digital balance (A & D Systems, Milpitas, 
California, USA) to obtain their fresh mass (FM). Leaf discs were then 
saturated to full turgor by immersing them in distilled water in the dark for 24 
h. Leaves were reweighed to obtain their turgid mass (TM) and put into an oven 
at 70°C for 3 d to obtain their dry mass (DM). RWC was calculated as: RWC = 
(FM — DM)(TM — DM) xX 100. 


Repeated physiological measurements—Carbon assimilation (A), stomatal 
conductance (g), intracellular CO2 concentrations (Ci), chlorophyll fluores- 
cence (F\/Fm), and spectral reflectance were measured on a subset of 40 
individuals (10 individuals from each species-treatment combination) at each of 
the sampling intervals (onset, 2 wk, 4 wk, 6 wk, and 7 wk) during the drought. 
A mature, fully expanded leaf was chosen at the onset of the experiment for the 
repeated physiological measurements, and the same leaf was used at each of the 
subsequent sampling intervals. All repeated physiological measurements were 
conducted between 0800 and 1130 hours. 

Maximum photosynthesis (Amax), g, and Ci were measured in situ using an 
open-system portable infrared gas analyzer (Welles, 1986) (LI-6400, Li-Cor, 
Lincoln, Nebraska, USA). Measurements of Amax were taken under saturating 
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P values for effects included in the analysis of variance (ANOVA) for the various physiological measurements of Quercus sideroxyla and 


Q. laceyi using mixed models for (A) repeated physiological traits and P values for one-way ANOVAs for (B) destructively sampled physiological 


traits. 


(A) Repeated measures physiological traits 


Root collar No. 

Fixed factor SIPI WBI chiNDI NDVI REP diameter Height leaves Amax g FV/Fin 
Treatment 0.0529 0.0008 0.0513 0.0001 0.0001 0.0017 0.0100 0.3632 0.0000 0.0000 0.0000 
Species 0.0141 0.8663 0.0002 0.0023 0.0040 0.2259 0.0260 0.0036 0.0000 0.0000 0.0012 
Time 0.0778 0.0001 0.0001 0.0001 0.0001 0.8178 0.8000 0.1058 0.0000 0.0000 0.0000 
(B) Destructively sampled physiological traits 

Shoot Root 

Fixed factor Poa aa RWC mass mass R:S 
Treatment 0.0001 0.0001 0.0120 0.0030 0.0001 0.8625 
Species 0.6489 0.0057 0.0010 0.5797 0.1263 0.3723 
Time 0.0001 0.0001 0.0001 n/a n/a n/a 
Treatment X time 0.0001 0.0200 0.0120 n/a n/a n/a 
Treatment X species 0.5385 0.5734 0.2000 0.8748 0.0843 0.4326 
Species X time 0.4726 0.0457 0.2530 n/a n/a n/a 
Species X time X treatment 0.9220 0.8220 0.1209 n/a n/a n/a 


Note: Amax = maximum photosynthesis; g = stomatal conductance; F\y/Fm = chlorophyll fluorescence to maximal chlorophyll fluorescence; chiNDI = 
chlorophyll normalized difference vegetation index; NDVI = normalized difference vegetation index; REP = red edge position; RWC = relative water 
content; SIPI = structure independent pigment index; WBI = water band index; ‘pa = pre-dawn water potential; ‘ma = mid-day water potential. 


light (1800 umol-m~?-s~!) using the LED attachment for the LI-6400 
photosynthesis system. Light curves were used to verify saturating light levels 
for the two species (data not shown). 

The maximum efficiency of photosystem II (PSII) chemistry (Fy/Fm) was 
measured in situ on attached leaves using a portable pulse-modulated 
fluorometer (OS 500, Opti-Sciences, Hudson, New Hampshire, USA). Leaves 
were dark-adapted for 1 h, which was sufficient to allow complete relaxation of 
energy-dependent quenching. The ratio F\/Fm was calculated as (Fm — Fo)/Fm, 
where Fim and F are the maximum and basal fluorescence yields, respectively, 
of dark-adapted leaves (Gentry et al., 1989). 

Spectral reflectance of the upper leaf surface (adaxial side) was measured 
using a UNIspec Spectral Analysis System (PP Systems, Haverhill, 
Massachusetts, USA) over the range of 400-1000 nm with a 2.0 mm diameter 
foreoptic and an internal 6.8 W halogen lamp. Standard controls were dark and 
a spectralon reflectance white standard. The UNIspec unit was recalibrated after 
every 10 measurements using the reflectance white standard. Leaves were held 
in a black polyvinyl chloride (PVC) clip at a 60° angle relative to the foreoptic 
for scanning (Rj, = leaf radiance at wavelength 4). We used four indices based 
on the reflectance properties of plant pigments and their biochemical 
components to identify differences in the reflectance spectra according to 
species and treatment: (1) a measure of carotenoids: chlorophyll a ratio, the 
structure independent pigment index (Pefiuelas et al., 1995) (SIPI) = (Rso0 — 
R4a4s)/(Rs00 — Rego); (2) the normalized difference vegetation index, NDVI = 
(R750 — Reso)/(R750 + Reso); (3) a revised version of NDVI that is sensitive to 
chlorophyll a concentrations (Gitelson and Merzlyak, 1996), chNDI = (R750 — 
R705)/(R750 + R705); and (4) a water band index, WBI = (Rooo + (Rooo — Rs30))/ 
Ro70, a measure of vegetation water content (Riedell and Blackmer, 1999). 

The red edge position (REP) is the name given to the change in reflectance 
in the 670-780 nm range of vegetation spectra caused by strong chlorophyll 
absorption by photosystems I and II and by internal leaf scattering. It is defined 
as the wavelength of the inflection point of the reflectance slope across this 
region and is indicative of leaf chlorophyll content. Under stress, the red edge 
position shifts toward the shorter wavelengths as a result of senescence or 
stress-induced chlorosis (Rock et al., 1988; Pefiuelas et al., 1994; Thorhaug 
et al., 2006). REP was calculated using the four-points formula developed by 
Guyot et al. (1992), which performs well in comparison to other more 
computationally intensive REP measures (Pu et al., 2003). The REP was 
determined by calculating the reflectance at the red edge inflection point (Rye) 
given by Rre = (Re70 + R780)/2, and then by deriving the red edge wavelength 
(Are) given by Are = 700 + 40*[(Rre — R700)/(R740 — R700)] where R670, R700, 
R740, and R7g0 are the reflectance values at 670, 700, 740, and 780 nm 
wavelength, respectively. 


Growth and biomass allocation—Seedling growth and dry mass gain were 
determined at the end of the experiment for 20 individuals in each species- 
treatment combination (N = 80) in a stratified random manner. Species height 
and leaf number were determined on the same individuals at the onset of the 
experiment and at wk 4 and 7. Root collar diameter was measured at wk 4 and 
7. Seedlings were harvested at the end of the experiment, and leaves, stems, and 
roots were separated, dried at 70°C for 3 d, and then dried at 90°C until a 
constant mass was reached. Whole shoot and root dry masses were recorded for 
each individual, and the root : shoot ratio (R :S) was calculated. The differences 
between root and shoot masses in the control and the drought treatment were 
also determined. 


Statistical analyses—Multifactor ANOVAs were used to determine the 
effects of species, treatment, and time on seedling performance. All statistical 
analyses were performed using the R Statistical Language (R Development 
Core Team, 2005). RWC, pa, ‘ma, dry plant matter mass (leaves, stems, 
roots), and root: shoot measurements were analyzed using a general linear 
model (GLM) procedure with fixed effects of species, treatment, and time and 
the interaction effects of species X treatment, species X time, treatment X time, 
and species X treatment X time. Differences between species—treatment 
combinations at each sampling interval were determined using a Tukey’s test of 
significant difference. Differences between ‘Ppa and ‘Wma measurements were 
also were assessed using a GLM. Repeated measurements of Amax, Ci, g, SIPI, 
PRI, ch,NDI, NDVI, WBI, RE, F\/Fm, and seedling growth (number of leaves, 
height, root collar diameter) were analyzed using separate ANOVAs with 
repeated measures by a mixed models procedure with species and treatment as 
fixed effects and time as a random effect. Mixed models were used rather than 
repeated measures GLMs to account for differences in the covariance structure 
of the repeated measures. A regression was also performed to assess the 
relationship between Amax and g. 


RESULTS 


Water relations—The ANOVAs (Table 1) revealed signif- 
icant (P < 0.0001) effects of treatment, time, and the 
interaction of treatment X< time for ‘Yma, which demonstrated 
the impact of the drought treatment on the seedlings of both 
species, and showed the accelerating impact of drought over 
time. Species-level differences for ‘Wma were significant at 
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TaBLe 2. Means by species and treatment for different (A) physiological measurements; (B) reflectance indices; and (C) growth measurements of 
Quercus sideroxyla (QUSI) and Q. laceyi (QULA). Means (+SE) are presented in relation to Table 1. Lowercase letters beside the SEs represent 


significant differences between species and treatments at the 5% level according to Tukey’s comparisons. 


(A) Physiological 


measurement Species Treatment Week 0 Week 2 Week 4 Week 6 Week 7 
QULA Control 670+ 19a 722+ 2.15a 67.1 + 1.lo6a 55.8 + 0.69 a 59.2 + 0.76 a 
QUSI Control 753+ 14b 72.5 + 0.97 a 76.5 + 1.63 b 66.5 + 1.10 a 63.9 + 0.13 a 
QULA Drought 618+ 1.7a 58.3 + 1.87 a 57.5 + 1.22a 44.0 + 0.96 a 44.0 + 2.78 a 
RWC (%) QUSI Drought 81.2 + 1.0b 74.8 + 1.50a 72.4 + 0.52 b 53:3:+°2:53-a 54.5 + 1.6la 
QULA Control —0.3125 + 0.01 a —0.37 + 0.01 a —0.28 + 0.02 a —0.2 + 0.00 a —1.31 + 0.21 ac 
QUSI Control —0.225 + 0.01 b —0.32 + 0.00 a —0.60 + 0.00 a —0.225 + 0.00 a —0.35 + 0.00 ac 
QULA Drought —0.35 + 0.03 a —0.46 + 0.02 a —0.60 + 0.00 b —2.85 + 0.25 b —3.35 + 0.35 ac 
pa (MPa) QUSI Drought —0.19 + 0.00 b —0.50 + 0.02 a —0.51 + 0.03 b —1.86 + 0.29 b —2.6 + 0.09 b 
QULA Control —0.516 + 0.13 a —0.88 + 0.38 a —0.70 + 0.00 a —0.6 + 0.17 a —1.7 + 1.99 ac 
QUSI Control 0.25 + 0.07 b —0.82 + 0.34 a —0.47 + 0.23 a —0.6 + 0.20 a —0.6 + 0.00 a 
QULA Drought —0.567 + 0.10 a —0.83 + 0.15 a -1.13 + 0.56 a —2.2 + 0.95 b —6.6 + 1.56 be 
Yina (MPa) QUSI Drought —0.542 + 0.04 b —0.86 + 0.05 a —1.00 + 0.95 a —-14 + 1.95b —4.5 + 0.76 b 
QULA Control 10.374 + 0.55 a 10.49 + 0.51 a 8.18 + 0.48 a 9.532 + 0.66 a 7.923 £0.71 a 
QUSI Control 15.610 + 0.48 b 14.29 + 0.57 b 10.57 + 0.70 a 13.442 + 0.63 a 11.817 + 0.56 a 
QULA Drought 4.083 + 0.33 a 4.96 + 0.31 a 4.20 + 0.66 b 1.194 + 0.19b 0.409 + 0.06 b 
Amax (jumol-m~?-s~!) QUSI Drought 14.324 + 0.76 b 12.39 + 0.53 b 8.59 + 0.60 b 1.847 + 0.47 b 0.758 + 0.33 b 
QULA Control 0.187 + 0.01 a 0.17 + 0.01 a 0.09 + 0.01 a 0.144 + 0.01 a 0.107 = 0.01 a 
QUSI Control 0.286 + 0.01 b 0.21 + 0.01 b 0.14 + 0.01 a 0.211 + 0.0la 0.180 = 0.01 a 
QULA Drought 0.244 + 0.01 a 0.12 + 0.0la 0.05 + 0.01 b 0.011 + 0.00 b 0.002 + 0.00 b 
g (umol-m-?-s~!) QUSI Drought 0.255 + 0.02 b 0.20 + 0.01 b 0.09 + 0.01 b 0.008 + 0.01 b 0.007 + 0.00 b 
QULA Control 0.838 + 0.00 a 0.83 + 0.00 a 0.82 + 0.00 a 0.822 + 0.00 a 0.782 + 0.02 a 
QUSI Control 0.850 + 0.00 a 0.85 + 0.00 b 0.83 + 0.01 a 0.810 + 0.01 b 0.850 + 0.00 a 
QULA Drought 0.842 + 0.00 a 0.81 + 0.0la 0.80 + 0.01 a 0.783 + 0.01 ac 0.276 + 0.05 b 
Fy/Fim QUSI Drought 0.832 + 0.01 a 0.84 + 0.00 b 0.81 + 0.01 a 0.525 + 0.06 be 0.265 + 0.06 b 
(B) Reflectance index Species Treatment Week 0 Week 2 Week 4 Week 6 Week 7 
QULA Control 0.776 += 0.00 a 0.799 + 0.00 a 0.799 + 0.00 a 0.743 = 0.02 a 
QUSI Control 0.827 + 0.01 b 0.809 + 0.02 a 0.809 + 0.02 a 0.721 + 0.04 a 
QULA Drought 0.780 + 0.00 a 0.792 + 0.00 a 0.792 + 0.00 a 0.685 + 0.02 a 
NDVI QUSI Drought 0.828 + 0.00 b 0.845 + 0.00 a 0.845 + 0.00 a 0.593 = 0.05 a 
QULA Control 0.452 + 0.00 a 0.438 + 0.01 a 0.438 + 0.01 a 0.424 + 0.01 a 
QUSI Control 0.543 + 0.01 b 0.527 + 0.02 b 0.527 + 0.02 b 0.536 + 0.00 a 
QULA Drought 0.445 + 0.00 a 0.454 + 0.00 a 0.454 + 0.00 a 0.310 + 0.02 b 
chiINDI QUSI Drought 0.500 + 0.01 b 0.509 + 0.01 b 0.509 + 0.01 b 0.260 + 0.03 b 
QULA Control 1.029 + 0.00 a 1.028 + 0.00 a 1.028 + 0.00 a 1.019 + 0.00 a 
QUSI Control 1.030 + 0.00 b 1.028 + 0.00 a 1.028 + 0.00 a 1.022 + 0.00 a 
QULA Drought 1.028 + 0.00 a 1.027 + 0.00 a 1.027 + 0.00 a 1.005 + 0.00 b 
WBI QUSI Drought 1.030 + 0.00 b 1.030 + 0.00 a 1.030 + 0.00 a 0.996 + 0.00 b 
QULA Control 1.033 + 0.00 a 1.028 + 0.00 a 1.028 + 0.00 a 0.994 + 0.00 ac 
QUSI Control 1.035 + 0.00 b 1.048 + 0.01 a 1.048 + 0.01 a 1.016 + 0.00 ac 
QULA Drought 1.030 + 0.00 a 1.025 + 0.00 a 1.025 + 0.00 a 1.005 + 0.01 b 
SIPI QUSI Drought 1.039 + 0.00 b 1.030 + 0.00 a 1.030 + 0.00 a 1.261 + 0.05 b 
QULA Control 715.7 + 0.20 a 715.6 + 0.21 a 716.1 + 0.16 a 715.3 = 0.20 a 
QUSI Control 718.4 + 0.19 b 717.6 = 0.52 b 718.3 + 0.26 a 717.8 £0.12 a 
QULA Drought 716.1 £0.15 a 716.2 + 0.18 a 716.3 + 0.19b 712.6 = 0.58 b 
REP (A) (nm) QUSI Drought T17.1 + 0.35 b 717.2 + 0.34b 715.2 + 0.47 b 712.0 = 0.93 b 


P < 0.01. Differences between ‘pa and ‘Wma were significant 
(P < 0.0001) according to a GLM test. Quercus sideroxyla had 
significantly higher Ppa, ma, and RWC than Q. laceyi 
(P < 0.01) at the onset of the experiment. Species differences 
in ‘Ypa and Wma continued throughout the drought, with Q. 
laceyi having consistently lower ‘Ppa and Wma than Q. 
sideroxyla (P_ < 0.05). Treatment-level differences in pa 
and Wma were significant (P < 0.05) beginning in wk 4 
(Table 2). 


Photosynthesis, stomatal conductance, intracellular CO2, 
and chlorophyll fluorescence—All ANOVAs were signifi- 
cant at P < 0.0001 (Table 1) for Amax, g, and F\/Fm. 
Intracellular CO2 (Ci) concentrations were not significantly 


different between species, treatments, or sampling intervals, 
and they had an overall mean of 112.1 (£7.95 SE) pmol 
CO2-mol air~!. Species differences in photosynthesis and 
stomatal conductivity were evident from the onset of the 
experiment (P < 0.001), and Amax and g were higher for Q. 
sideroxyla than for Q. laceyi throughout the experiment 
(Table 2). Species-level differences in Amax and g declined in 
latter weeks of the drought, while treatment-level differences 
increased, with significant differences (P < 0.05) occurring 
from the 4th wk on. Amax values were lower for Q. laceyi 
drought treatment seedlings during the first two sampling 
intervals due to random variability. This was probably the 
result of this species’ existence on a wider range of sites 
relative to Q. sideroxyla. A nonsignificant decline in Amax 
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TABLE 2. Continued. 
(C) Growth measure Species Treatment Week 0 Week 2 Week 4 Week 6 Week 7 
QULA Control 17.67 + 0.87 a 22.22 + 0.94 a 18.44 + 1.03 a 
QUSI Control 20.30 + 1.01 a 23.50 + 1.49 a 21.30 + 0.77 a 
QULA Drought 16.27 + 0.98 a 18.09 + 1.22b 15.09 + 0.66 a 
No. leaves QUSI Drought 16.80 + 0.62 a 18.30 + 1.17b 17.42 + 1.40a 
QULA Control 17.67 + 0.94 a 22.22 + 0.94 a 18.44 + 1.03 a 
QUSI Control 20.30 + 1.49 a 23.50 + 1.49 a 21.30 + 0.77 a 
QULA Drought 16.27 + 1.22 a 18.09 + 1.22b 15.09 + 0.66 b 
Height (cm) QUSI Drought 16.80 + 1.17 a 18.30 + 1.17b 17.42 + 1.40b 
QULA Control 4.95 + 0.56 a 3.54 + 1.95 b 
QUSI Control 4.37 + 0.02 a 3.16 + 0.56 a 
QULA Drought 3.42 + 0.45 a 3.68 + 0.46 a 
Root collar (cm) QUSI Drought 2.88 + 0.02 a 4.70 + 0.23 b 
QULA Control 144+ 0.19a 
QUSI Control 1.93 + 0.19 a 
QULA Drought 1.09 + 0.08 b 
Root mass (g) QUSI Drought 1.06 + 0.08 b 
QULA Control 0.43 + 0.08 a 
QUSI Control 0.45 + 0.04 a 
QULA Drought 0.25 + 0.03 b 
Shoot mass (g) QUSI Drought 0.29 + 0.03 b 
QULA Control 4.74 + 0.57 a 
QUSI Control 4.69 + 0.50 a 
QULA Drought 5.25 + 0.56a 
Root : shoot QUSI Drought 4.36 + 0.45 a 


Note: Amax = maximum photosynthesis; g = stomatal conductance; F\/Fm = chlorophyll fluorescence to maximal chlorophyll fluorescence; chiINDI = 
chlorophyll normalized difference index; NDVI= normalized difference vegetation index; REP = red edge position; RWC = relative water content; SIPI = 
structure independent pigment index; WBI = water band index; ‘¥pa = pre-dawn water potential; ‘Yma = mid-day water potential. 


occurred in control seedlings beginning on wk 4; the decline 
was probably the result of leaf aging during the experiment 
because all repeated photosynthesis measurements were taken 
on the same leaf. A positive relationship existed between Amax 
and g for the pooled data set (Fig. 2), indicating an increase in 
stomatal conductance with increased photosynthesis. The 
ratios of Fy/Fm were much less sensitive than Amax and ‘Y, and 
they differed little between treatments or species until wk 6 
and 7 of the drought (Table 2). Chlorophyll fluorescence in Q. 
sideroxyla declined first during wk 6 of the drought, while Q. 
laceyi had no significant decline in chlorophyll fluorescence 
until wk 7 (P < 0.05). 


Spectral reflectance—Reflectance spectra varied between 
species, although spectra for both species were characterized 
by a broad peak at 550 nm, a trough at 680 nm, a sharp increase 
through the red edge around 700 nm, followed by another 
broad peak at 760 nm, and a gently decreasing plateau above 
750 nm (Fig. 3). Averaged across all control treatment 
seedlings and time intervals, Q. laceyi reflectance was higher 
than Q. sideroxyla for all wavelengths. Differences between 
control and drought treatments were observed 6 and 7 wk after 
drought onset (Fig. 3). The drought resulted in increased 
reflectance for both species relative to control treatments at all 
wavelengths and in a broad peak in reflectance that ranged 
across the visible light spectrum (400-700 nm). 

Reflectance ratios were calculated to illustrate the differences 
in the reflectance spectra between species and treatments for wk 
6 and 7 of the drought (Fig. 4). The reflectance ratios were 
calculated using the mean reflectance of each species—treatment 
combination, and they included drought Q. laceyi: control Q. 
laceyi (DL: CL); drought Q. sideroxyla: control Q. sideroxyla 
(DS: CS); drought Q. sideroxyla: drought Q. laceyi (DS : DL); 


control Q. sideroxyla:control Q. laceyi (CS:CL). Major 
differences were similar for wk 6 and 7 of the drought. The 
two species’ drought to control ratios were the highest (DS : CS 
and DL: CL) with peaks at 660 and 480 nm, and the change in 
reflectance due to drought was greater for Q. sideroxyla than for 
Q. laceyi. Species differences in reflectance in drought and 
control treatments were similar (DS:DL and CS:CL), with 
increases in reflectance near 680 nm. 


Reflectance indices—The treatment effect was significant 
(P < 0.05) for all reflectance indices (Table 1). Species effects 
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Fig. 2. Relationship between stomatal conductance (g) and maximum 
photosynthesis (Amax) for pooled samples. 
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Reflectance spectra of species—treatment combinations 6 wk (A) and 7 wk (B) after the onset of the drought experiment. Abbreviations: CL, 


control, Quercus laceyi; CS, control, Q. sideroxyla; DL, drought, Q. laceyi; DS, drought, Q. sideroxyla. 


were significant (P < 0.05) for SIPI, cnNDI, and NDVI but not 
for WBI. Time since watering was also significant for all 
indices except SIPI. Species, treatment, and time effects were 
evident in chiNDI values (P < 0.05) (Table 2). Like the other 
measures of seedling performance, the reflectance indices 
identified significant differences between species at the 
beginning of the experiment, with Q. sideroxyla having 
significantly (P < 0.05) higher values than Q. /aceyi for all 
indices (Table 2). Species-level differences diminished over the 
course of the experiment and were replaced by treatment-level 
differences (P < 0.05). 

The red edge A differed significantly (P < 0.01) by species, 
treatment, and time (Tables 1 and 2). Quercus sideroxyla had a 
higher REP than Q. /aceyi at all samplings until wk 6 of the 
drought, when treatment-level differences became significant 
(P < 0.05). The shift in REP toward shorter wavelengths in 
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reflectance ratio 
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drought seedlings began in wk 6, although treatment effects 
were most pronounced in wk 7. 


Growth and biomass allocation—Treatment effects on plant 
height were evident in wk 4 and 7 (P < 0.05) (Tables 1 and 2). 
Quercus sideroxyla had significantly more leaves than Q. 
laceyi (Tables 1 and 2), but there were no significant treatment 
effects on leaf number. Control seedlings had significantly 
(P < 0.001) thicker root collar diameters than drought 
seedlings (Tables 1 and 2). 

Control seedlings had significantly greater (P < 0.01) root 
and shoot biomass than drought seedlings, and control Q. 
sideroxyla had greater (P < 0.001) dry root mass than Q. laceyi 
(Table 1). Decline in root and shoot mass was significantly 
greater (P < 0.01) with Q. sideroxyla than Q. laceyi (Fig. 5A, 
B). Quercus sideroxyla roots were longer and more branched 
than those of Q. laceyi, and control seedling roots were longer 


week 7 
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Fig. 4. Reflectance ratio spectra for species—treatment combinations 6 wk (A) and 7 wk (B) after the onset of the drought experiment. Abbreviations: 
CL, control, Quercus laceyi; CS, control, Q. sideroxyla; DL, drought, Q. laceyi; DS, drought, Q. sideroxyla. 
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Mean differences (+SE) between control and drought shoot (A) and root (B) biomass of seedlings of Quercus laceyi and Q. sideroxyla at the 


end of the 7-wk drought. Control and drought treatments differed significantly according to Tukey’s comparisons (P < 0.01). 


and more branched than drought seedling roots. Root: shoot 
ratios did not differ significantly between species or treatments. 


DISCUSSION 


Our results suggest that the oaks of northern Mexico have 
evolved to survive and reproduce under a range of climatic 
conditions. The species-level differences in oak physiology 
reflect each species’ adaptation to local growing conditions, 
potentially explaining the differentiation and diversification of 
oaks within and between mountain ranges in Mexico. The 
similar drought responses of Q. laceyi and Q. sideroxyla 
suggest that both species are well adapted to drought, although 
their performance differs substantially under conditions that are 
not moisture-limited. 

The dominance of Q. sideroxyla on cooler, wetter sites is 
explained by its higher photosynthesis, stomatal conductance, 
relative water content, and water potentials and its lower 
spectral reflectance relative to Q. laceyi, conferring competitive 
advantages on high-productivity sites at upper elevations. 
Likewise, the dominance of Q. laceyi on hotter, drier sites 
reflects its adaptation to the more extreme conditions typical 
at lower elevations. Our water potential (Vpa and ma) and 
reflectance ratio data corroborate these findings, by showing 
that Q. laceyi was better able to carry out key physiological 
processes than was Q. sideroxyla under drought conditions. 
This suggests that Q. sideroxyla avoids drought as a high- 
elevation site specialist, while Q. laceyi is a drought tolerant 
species that competes well under more extreme, low elevation 
conditions. 

In general, higher water potential, higher photosynthetic rate, 
and higher relative water content equate to less drought 
tolerance in oaks (Nardini et al., 1999; Fotelli et al., 2000; 
Thomas and Gausling, 2000). While our results suggest that Q. 
sideroxyla is less drought adapted than is Q. laceyi, Q. 
sideroxyla remained highly resilient to drought, even with its 
initially higher water potentials, photosynthetic rate, and 
relative water content. A possible explanation for this pattern 


is that the high variability in year-to-year precipitation in 
southwestern North America has caused species across a range 
of elevations to be adapted to wet and dry climatic cycles in 
order to persist long-term. Such findings further suggest that 
even less drought tolerant species such as Q. sideroxyla may be 
able to persist under future climate change scenarios (Hough- 
ton et al., 2001). 


Biomass allocation and growth—The first and most 
sensitive response to water deficiency is a decrease in turgor 
and a slowing of the growth processes (Larcher, 2003). Plants 
of both species were shorter, had smaller root collars, and 
allocated less biomass to roots and shoots in the drought than in 
the control treatment because the rate of cell division and the 
duration of leaf expansion drop as cell turgor declines. 
Likewise, the drop in photosynthesis during the drought 
limited COz assimilation, thereby limiting growth. 


Water relations—Plant water status is controlled by the 
relative rates of water absorption (uptake by roots) and water 
loss (transpiration), and a plant’s drought tolerance depends 
upon its ability to tolerate internal water deficits. As soil 
moisture is depleted under drought, the leaf water content and 
the overnight recovery of are inhibited by progressive 
increases in liquid flow resistance by the soil or plant 
(Kozlowski and Pallardy, 1997). The lower RWC and water 
potentials (Ppa and Wma) of Q. laceyi relative to Q. sideroxyla 
suggest that Q. laceyi is more drought tolerant than Q. 
sideroxyla because it is better able to endure extreme water 
deficits, although these differences in adaptation did not 
translate to species’ differences in growth or biomass. 

At wk 7 in the drought treatment, ‘Wpa and Wma of both 
species were more negative (Q. laceyi = —6.6 MPa and Q. 
sideroxyla = —4.5 MPa) than water potentials of seedlings of 
oak species from other geographical regions including Q. suber 
(—3.0 MPa) in Spain (Aranda et al., 2005) and Italy (-1.5 MPa) 
(Nardini et al., 1999), QO. robur (—3.0 MPa) and Q. petraea 
(—3.2 MPa) in Germany (Thomas and Gausling, 2000), Q. 
rubra (—2.5 MPa) in the northeastern United States (Kubiske 


816 AMERICAN JOURNAL OF BOTANY 


and Abrams, 1992; Cavendar-Bares, 2000), and several other 
oak species of the southwestern United States (Q. buckleyi, Q. 
lobata, Q. macrocarpa, Q. muehlenbergii, and Q. shumardii) 
(Balok and St. Hilaire, 2002). 

Week 7 ‘Ppa and ‘Yma compare more closely to other oaks, 
pinyon pines, and junipers from other arid regions of the 
Southwest, and oaks such as Q. gambelii and Q. turbinella 
appear to tolerate even lower ‘Wma (—7.3 MPa and —7.0 MPa, 
respectively) upon exposure to a similar 7-wk drought than did 
the seedlings in our study (Neilson and Wullstein, 1987). This 
suggests that oak trees from the southwestern United States and 
northern Mexico are more drought tolerant than oaks from 
other geographical regions, but that there is a lower limit of 
~—7.0 MPa for oak water stress tolerance. 


Photosynthesis—The decreases in ‘Ppa and ‘Pima in wk 6 and 
7 of the drought induced stomatal closure, which decreased the 
supply of CO2 to the mesophyll cells, and consequently, 
decreased the rate of leaf photosynthesis. Increases in internal 
water deficit are associated with disruption of biosynthetic 
activities, including the generation of the chlorophyll necessary 
for photosynthesis (Fitter and Hay, 2002), and photophosphor- 
ylation and photosynthetic carbon metabolism may be 
particularly sensitive to dehydration because of the toxic 
effects of Mg?* concentrations that accompany the removal of 
water from chloroplasts in dehydrating leaves (Kozlowski and 
Pallardy, 1997). 

The higher photosynthetic rates of Q. sideroxyla relative to 
Q. laceyi also correspond to their elevational distributions. 
Adiabatic lapse rates and increases in relative humidity with 
increasing elevation can ameliorate climates in xeric environ- 
ments (Barry, 1992; Robinson et al., 1999; Whiteman, 2000), 
providing more favorable growing conditions for plants at 
higher elevations. The pattern of increased photosynthesis with 
increasing elevation in our study is the reverse of the patterns 
observed in other alpine environments near upper treeline 
(Tranquillini, 1979). In the Southwest, the harshest sites are at 
low elevations, where lower treelines are limited by moisture 
and temperature (VanDevender, 1995). Lower temperatures 
and higher relative humidity at high elevations can therefore 
explain the pattern of high photosynthesis by Q. sideroxyla 
relative to Q. laceyi, and the lower maximal photosynthesis of 
Q. laceyi is likewise explained by its distribution in a more 
water-limited environment that demands lower maximum 
metabolic rates. 


Chlorophyll fluorescence—tThe reduction in F\y/Fm during 
the drought reflects a decline in the efficiency of light use for 
electron transport by PSII. Under drought conditions, the 
decreases in photosynthesis are accompanied by a decrease in 
the rate of consumption of ATP and NADPH for COz 
assimilation, which can reduce the rate of linear electron 
transport from the quinone acceptors of PSII, and conse- 
quently, in F\y/Fm (Baker and Rosenqvist, 2004). The decline in 
F\/Fm of Q. sideroxyla in wk 6 of the drought and a similar 
decline by Q. laceyi in wk 7, support our findings from the 
water relations data that low elevation Q. laceyi is more 
drought tolerant than high elevation Q. sideroxyla. These data 
suggest that Q. laceyi can better compensate for decreases in 
the use of ATP and NADPH in photosynthesis metabolism 
resulting from drought. 
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Spectral reflectance—Leaf optical properties are well 
correlated with vegetation water content (Sims and Gamon, 
2003) and leaf pigment content (Curran, 1989; Gamon and 
Surfus, 1999; Sims and Gamon, 2002), and reflectance is 
known to shift in response to environmental stress (Pefiuelas 
et al., 1994; Rock et al., 1998; Carter and Knapp, 2001; 
Thorhaug et al., 2006). The higher reflectance spectrum of Q. 
sideroxyla, along with the higher chiNDI index and reflectance 
ratios of control seedlings, suggest that Q. sideroxyla has a 
higher concentration of plant pigments than Q. laceyi, which 
relates to its higher maximum photosynthetic capacity. 

Changes in the reflectance properties of drought-stressed 
seedlings are the product of decreases in plant water balance 
and photosynthetic pigment concentrations. The drop in WBI 
in latter weeks of the drought indicated decreased plant RWC 
(Pefiuelas et al., 1997; Sims and Gamon, 2003). Decreases in 
chINDI, increases in SIPI, the shift of the REP to shorter 
wavelengths, and the increases in reflectance observed over a 
broader portion of the visible spectrum are typical plant 
responses to stress (Pefiuelas et al., 1994; Carter and Knapp, 
2001; Richardson et al., 2001; Richardson and Berlyn, 2002; 
Thorhaug et al. 2006), although the position of the red edge 
was high compared to other studies (i.e., Richardson et al., 
2001; Richardson and Berlyn, 2002). Furthermore, the 
reflectance ratios between species in the drought and control 
reflectance spectra had noticeable peaks in the 680 nm region, 
suggesting a greater shutdown of PSII in Q. sideroxyla than in 
Q. laceyi during wk 6 and 7. 

While photosynthetic rates tended to increase with oak 
elevational distribution, the reflectance spectra and reflectance 
ratios show an opposite trend, with lowland Q. laceyi having 
higher reflectance than upland Q. sideroxyla. This pattern is 
consistent with work in mountains of the northeastern United 
States (Richardson et al., 2001; Richardson and Berlyn, 2002), 
and the observed differences in spectral reflectance also 
probably relate to variation in the two species’ growing 
conditions. The reflectance data suggest that the light green 
leaves of Q. laceyi have lower chlorophyll concentrations, 
which correspond both to its lower photosynthetic rate and its 
distribution on hot, dry, low-elevation sites. In contrast, the 
lower reflectance spectrum of Q. sideroxyla corresponds to 
higher chlorophyll absorption and higher photosynthetic 
capacity, most likely related to its distribution on cooler, 
wetter, high-elevation sites. 


Application of greenhouse results to field conditions—The 
application of greenhouse experiments to field conditions 
warrants discussion (e.g., Harper, 1977), because the field 
environment can differ dramatically from controlled green- 
house conditions. We addressed this by mimicking natural 
patterns of oak germination, growth, and drought timing in the 
southwestern United States and northern Mexico. 

Results from this study correspond to other field-based 
studies of plant water potential in the southwestern United 
States (i.e., Lajtha and Barnes, 1991; Lajtha and Getz, 1993; 
Linton et al., 1998), suggesting that our experiment approx- 
imated field soil moisture conditions. Seedling survival over 
the course of the acute drought was similar to other acute 
drought studies (Neilson and Wullstein, 1987; Lajtha and 
Barnes, 1991; Barton and Teeri, 1993; Lajtha and Getz, 1993; 
Linton et al., 1998; Balok and St. Hilaire, 2002). Although pots 
are known to constrain root elongation and subsequent water 
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uptake, roots of the seedlings in this study failed to reach the 
bottom of their pots during the experiment. 


Conclusions 


Results from this study demonstrate the physiological 
mechanisms responsible for the observed elevational separation 
and niche partitioning of Q. sideroxyla and Q. laceyi in the 
Sierra del Carmen, Coahuila, Mexico. Quercus laceyi is a 
drought tolerant species that can grow under the high water 
stress conditions typically observed at low elevations during 
summer months in northern Mexico, while Q. sideroxyla is a 
drought avoider species that grows poorly under extreme water 
stress and is therefore restricted to higher elevations with lower 
temperatures and higher relative humidity. However, both 
species are adapted to dry conditions, suggesting that they can 
survive under a range of climate scenarios. The differentiation 
of these two oaks according to elevation suggests that similar 
mechanisms may be responsible for the diversification of the 
genus Quercus across Mexico and the rest of North America. 
Future studies that investigate the niche partitioning of oaks 
across other mountain ranges in Mexico are warranted to 
explain the high diversity of this genus across Mexico. 
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